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Abstract 

We address the general question of how to reconstruct the field content 
of a quantum field theory from a given scattering theory in the context of 
the form factor program. For the 5C/(3)2-homogeneous Sine-Gordon model 
we construct systematically all n-particle form factors for a huge class of 
operators in terms of general determinant formulae. We investigate how 
different operators are interrelated by the momentum space cluster property. 
Finally we compute several two-point correlation functions and carry out the 
ultraviolet limit in order to identify each operator with its corresponding 
partner in the underlying conformal field theory. 
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1 Introduction 



The central concepts of relativistic quantum field theory, like Einstein causality 
and Poincare covariance, are captured in local field equations and commutation 
relations. As a matter of fact local quantum physics (algebraic quantum field the- 
ory) JIJ takes the collection of all operators localized in a particular region, which 
generate a von Neumann algebra, as its very starting point (for recent reviews see 
e.g. §). 

On the other hand, in the formulation of a quantum field theory one may alter- 
natively start from a particle picture and investigate the corresponding scattering 
theories. In particular for 1+1 dimensional integrable quantum field theories this 
latter approach has proved to be impressively successful. As its most powerful 
tool one exploits here the bootstrap principle || , which allows to write down exact 
scattering matrices. Ignoring subtleties of non-asymptotic states, it is essentially 
possible to obtain the latter picture from the former by means of the LSZ-reduction 
formalism ||]. However, the question of how to reconstruct at least part of the field 
content from the scattering theory is in general still an outstanding issue. Recently 
a link between scattering theory and local interacting fields in terms of polarization- 
free generators has been developed ||. Unfortunately, they involve subtle domain 
properties and are therefore objects which concretely can only be handled with 
great difficulties. 

In the context of 1+1 dimensional integrable quantum field theories the iden- 
tification of the operators is based on the assumption, dating back to the initial 
papers ||, that each solution to the form factor consistency equations || |7], f| |J 
corresponds to a particular local operator. Based on this numerous authors ||, |7|, ||, 
||, [T(| [ll], [H, [H, [14], ^] have used various ways to identify and constrain the specific 
nature of the operator, e.g. by looking at asymptotic behaviours, performing per- 
turbation theory, taking symmetries into account, formulating quantum equations 
of motion, etc. Our analysis will make especially exploit the conjecture that each 
local operator has a counterpart in the ultraviolet conformal field theory. 

In the present manuscript we show for a concrete model, the SU(3)2-homogeneous 
Sine-Gordon model (HSG), that, by means of the form factor program, it is possible 
to reconstruct the field content starting from its scattering matrix. Our analysis is 
based on the assumption [^, [U| that each solution to the form factor consistency 
equations || 0, f| |] corresponds to a particular local operator. We take furthermore 
into account that the SU(3)2-HSG model, like numerous other 1+1 dimensional in- 
tegrable models, may be viewed as a perturbed conformal field theory whose entire 
field content is well classified. Assuming now a one-to-one correspondence between 
operators in the conformal and in the perturbed theory, we can carry out an identifi- 
cation on this level, that is we associate each solution of the form factor consistency 
equations a local operator which is labeled according to the the ultra-violet confor- 
mal field theory. We therefore construct systematically all possible solutions for the 
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n-particle form factors related to a huge class of operators in terms of some general 
building blocks which consist out of determinants of matrices whose entries are ele- 
mentary symmetric polynomials depending on the rapidities. We demonstrate how 
these general solutions are interrelated by the momentum space cluster property. In 
particular we show that the cluster property serves also as a construction principle, 
in the sense that from one solution to the consistency equations we may obtain a 
huge class, almost all, of new solutions. Finally we compute the corresponding two- 
point correlation functions and carry out the ultraviolet limit in order to identify 
the corresponding conformal dimensions. 

Our manuscript is organized as follows: In section 2 we recall [HJ the solutions 
for the minimal form factors and the recursive equation which is central for the de- 
termination of the form factors. We describe the general structure of the n-particle 
form factors for a huge class of operators. In section 3 we provide a rigorous proof 
for all solutions. In section 4 we investigate the cluster property. In section 5 we 
compute several two-point correlation functions and carry out the ultraviolet limit 
on the base of a sum rule and the explicit two-point correlation function in order 
to identify the conformal dimensions of each operator. We state our conclusions 
in section 6. The appendix contains a collection of useful properties of elementary 
symmetric polynomials and some explicit formulae for the first non- vanishing form 
factors. 



2 The SU(3) 2 -HSG model form factors 

The SU(3)2-HSG model contains only two self-conjugate solitons which we denote, 
following the conventions of fllB] , by "+" and "— ". The two particle scattering 



matrix as a function of the rapidity 9 related to this model was found [17j to be 
S±± = —1 and S± T (9) = ± tanh - I 9 ± a — i— 

Here a is a real constant and corresponds to a resonance parameter. The system 
([1]) constitutes probably the simplest example of a massive quantum field theory 
involving two particles of distinct type. Nonetheless, despite the simplicity of the 
scattering matrix we expect to find a relatively involved operator content, since for 
finite resonance parameter the SU(3)2-HSG model describes a WZNW-coset model 
with central charge c = 6/5 perturbed by an operator with conformal dimension 
A = 3/5. Since this is true as long as a is finite, we shall be content in the 
following mostly by setting a = 0. It is expected from the classical analysis that for 
finite value of a we always have the same ultraviolet central charge and therefore the 



same operator content. The TBA-analysis carried out in |18] supports this analysis. 
Thus, a finite variation of a at the ultraviolet fixed point and away from it is not 
very illuminating and we therefore only distinguish the behaviour a — > oo and o 
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finite. In the former case one trivially observes that the S-matrices S±^ tend to 
one and the theory decouples into two Ising models. The related form factors have 
to respect this behaviour and all combinations involving different types of particles 
vanish, see [|TJ|. One may see also directly that the form factor solutions behave 
this way by employing the Riemann-Lebesgue theorem. 

The underlying conformal field theory has recently |19 found an interesting 



application in the context of the construction of quantum Hall states which carry 
a spin and fractional charges. 

Taking the scattering matrix as an input, it is in principle possible to compute 
form factors, by solving certain consistency equations || 0, f§ [|, an< ^ thereafter to 
evaluate correlation functions. Form factors are tensor valued functions, represent- 
ing matrix elements of some local operator 0{x) located at the origin between a 
multiparticle in-state and the vacuum, which we denote by 

F^~^{0 U . . . A) := (010(0)1^(00^(0,,) • • • y^n)). n ■ (2) 

Here the V^{6) are some vertex operators representing a particle of species /i. We 
commence now by recalling the basic ansatz for solutions of the form factors for the 



SU(3) 2 -HSG model from pq . We used the parameterization 

Z X ± m X =F 



Fn lfil ' ' ' "^(0! ...d n ) = H^"*»Q^"«»(x x . ...r,,)!! /••"/' (^,) (3) 



i<j 



where F^^ (%) := F^ 3 (6ij)/ (xf l + ^ 3 ■ The rapidities enter through the 

variable Xi = exp(0j) and the functions F^ 3 (0y) denote the so-called minimal form 
factors. They were found to be 

F£(B) = -i sinh | (4) 
m0) = 2^^-fexp(-/ f^F) ) = AF£ { Q. (5, 
Here G is the Catalan constant. For the overall constants we obtained 

j_^0\2s+T,m jsmi^s^s— m— 1+2t) ^srrw/2 jj-0\T,m ^_ Q 1 ^5) 

Note that at this point an unknown constant, that is H°' T ' m , enters into the proce- 
dure. This quantity is not constrained by the form factor consistency equations and 
has to be obtained from elsewhere. The polynomials Q have to satisfy the recursive 
equations 

Q ^ ( X, . . . , X n ) Dfl (x, Xi, . . . , X n )Q ^ (x'i ; • • • j 2>n) (7) 

-r m 

D l r(x, xx, ... , x n ) = k-ix) l+1 a+ £ x~ h (l - (-l)' +fc+ > fe • (8) 

/ fc=0 
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In particular 



D 



2s+T,2t+r' 



Ou • 30 - 



Xr, 



-i) 2s+T+1 a 



2sH 



-T /_ 4 

p=0 



2p+C 



(9) 



Here $ is related to the factor of local commutativity to = (— 1)^ = ±1. We 
introduced also the function ( which is or 1 for the sum d + t being odd or even, 
respectively. We shall use various notations for elementary symmetric polynomials 
(see appendix for some essential properties). We employ the symbol when the 
polynomials depend on the variables Xi, the symbol when they depend on the 
inverse variables x^ 1 , the symbol &k when they depend on the variables Xie~ a+t7T / 2 
and (Tfc when we set the first two variables to x\ = —x, X2 = x. The number of 
variables the polynomials depend upon is defined always in an unambiguous way 
through the l.h.s. of our equations, where we assume the first I variables to be 
associated with /x = + and the last m variables with // = — . In case no superscript 
is attached to the symbol the polynomials depend on all m + I variables, in case 
of a "+" they depend on the first I variables and in case of a "— " on the last m 
variables. 

Solving recursive equations of the type (0) in complete generality is still an 
entirely open problem. Ideally one would like to reach a situation similar to the one 
in the boostrap construction procedure of the scattering matrices, where one can 
state general building blocks, e.g. particular combinations of hyperbolic functions 
whenever backscattering is absent ^(J , infinite products of gamma functions when 
backscattering occurs or elliptic functions when infinite resonances are present. At 
least for all operators in the model we consider here this goal has been achieved. It 
will turn out that all solutions to the recursive equations (^) may be constructed 
from some general building blocks consisting out of determinants of matrices whose 
entries are elementary symmetric polynomials in some particular set of variables. 
Let us therefore define the (t + s) x (t + s)-matrix 



a 



2(J-i)+n 



i.l 



2{j-i)+2t+v 

More explicitly the matrix A reads 



AW 



for 
for 









1-4 




h6 
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1-2 
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h4 




















K+2 




1-4 




1-6 
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<v 


1-2 


<v 


h4 



















1 < i < t 
t <i < s + t 



\ 




(10) 



a 



2s+[i 






a 



2t+u / 
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The superscripts //, v may take the values and 1 and the subscripts /, m charac- 
terize the number of different variables related to the particle species "+", "— ", 
respectively. The different combinations of the integers (A,v,l,m will correspond to 
different kind of local operators O. 

In addition, the form factors will involve a function depending on two further indices 
Jl and v 

T, i~> , \ l — fn+U v — m , . 

C : =(^)~KrJ^ • (12) 

Here the /2, v are integers whose range, unlike the one for /i, u, is in principle not 
restricted. However, it will turn out that due to the existence of certain constraining 
relations, to be specified in detail below, it is sufficient to characterize a particular 
operator by the four integers fj,,u,l,m only. Then, as we shall demonstrate, all 
Q-polynomials acquire the general form 

We used here already a parameterization for l,m which will turn out to be most 
convenient. The subscripts in g and A are only needed in formal considerations, but 
in most cases the number of particles of species "+" and "— " are unambiguously 
defined through the l.h.s. of our equations. This is in the same spirit in which 
we refer to the number of variables in the elementary symmetric polynomials. We 
will therefore drop them in these cases, which leads to simpler, but still precise, 
notations. To illustrate this with examples, we consider for instance the solutions 
to the recursive equations ([?[) related to the trace of the energy momentum tensor 
0, the order operator E and the disorder operator /i which were already stated in 



g0|2 S +2,2 t+2 = ^Me-ft+DVj^^deti 1 - 1 (14) 

gE|2 S ,24+l = . S (2 t +3) g-hl detA 0,l ^ 
Q ^s,2t = .2.(4+1) ^-1,1 det ^0,0 _ (16) 

Here A is always taken to be a (t + s) x (t + s)-matrix. Notice that in comparison 
with (|T3| ) the factor of proportionality in ( |I~5|) and (Tl6| ) is only a constant, whereas in 
([y]) also the term <J\0\ appears. Terms of this type may always be added since they 
satisfy the consistency equations trivially. This is also the reason why in comparison 
with |16| we can safely drop in Q E l 2s > 2m the factor (a±) 2 (a^f) _ 2 . Additional reasons 



for this modification will be provided below. For O we were forced |T(| to introduce 
the factor 0\0\ in order to recover the solution of the thermally perturbed Ising 
model for 2s + 2 = 0. Note that for the value s — — 1 formally makes sense. 



3 Solution procedure 

We shall now recall the principle steps of the general solution procedure for the 
form factor consistency equations || [f[ ||, |§. For any local operator O one may 
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anticipate the pole structure of the form factors and extract it explicitly in form of 
an ansatz of the type @. This might turn out to be a relatively involved matter 
due to the occurrence of higher order poles in some integrable theories, e.g. [O], but 
nonetheless it is possible. Thereafter the task of finding solutions may be reduced to 
the evaluation of the minimal form factors and to solving a (or two if bound states 
may be formed in the model) recursive equation of the type (0). The first task can 
be carried out relatively easily, especially if the related scattering matrix is given as 
a particular integral representation ||. Then an integral representation of the type 
(H) can be deduced immediately. The second task is rather more complicated and 
the heart of the whole problem. Having a seed for the recursive equation, that is the 
lowest non-vanishing form factorf], one can in general compute from them several 
form factors which involve more particles. However, the equations become relatively 
involved after several steps. Aiming at the solution for all n-particle form factors, it 
is therefore highly desirable to unravel a more generic structure which enables one 



to formulate rigorous proofs. Several examples pi] , |], |ll[] have shown that often 
the general solution may be cast into the form of determinants whose entries are 
elementary symmetric polynomials. Presuming such a structure which, at present, 
may be obtained by extrapolating from lower particle solutions to higher ones or 
by some inspired guess, one can rigorously formulate proofs as we now demonstrate 
for the SU(3)2-HSG-model, for which some solutions were merely stated in |16| . 



We have two universal structures^ at our disposal. We could either exploit 
the integral representation for the determinant A, as presented in |16||, or exploit 
simple properties of determinants. Here we shall pursue the latter possibility. For 
this purpose it is convenient to define the operator Cfj (Rfj) which acts on the j th 
column (row) of an (n x n)-matrix A by adding x times the 2 th column (row) to it 

CfjA : A kj ^ A kj + xA hi l<i,j,k<n (17) 
Ri tj A : Ajk i-> A jk + xAik 1 <i,j,k <n . (18) 

Naturally the determinant of A is left invariant under the actions of Cfj and Rfj 
on A, such that we can use them to bring A into a suitable form for our purposes. 
Furthermore, it is convenient to define the ordered products, i.e. operators related 
to the lowest entry act first, 

C x a , b := II C; >p+1 and := {[ R^ . (19) 

p=a p=a 

*For the case at hand this is provided for some operators by the well known solutions of the 
Ising model. In general this is also a difficult hurdle to take as, for instance, one might need to 
know vacuum expectation values. 

t There exist also different types of universal expressions like for instance the integral repre- 
sentations presented in |J. However, these type of expressions are sometimes only of a very 
formal nature since to evaluate them concretely for higher n-particle form factors requires still a 
considerable amount of computational effort. 
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It will be our strategy to use these operators in such a way that we produce as 
many zeros as possible in one column or row of a matrix of interest to us. In order 
to satisfy (0) we have to set now the first variables in A to x\ = —x, %2 = x, which 
we denote as A thereafter and relate the matrices *4f+2,m an d ^-im- Taking relation 
fl70|) for the elementary symmetric polynomials into account, we can bring Af+ 2m 
into the form 



^n+2,s+t+i^i,s+t-iA 



l+2,m 



°2U-^ l<*<t 

a 2(j-i)+2t+u t <i < S + t ^ 

£ x^a- {p _ s _ 1)+u t = s + t + l 
{ p=i 



It is now crucial to note that since the number of variables has been reduced by two, 
several elementary polynomials may vanish. As a consequence, for 2s + 2 + > I 
and 2t + 2 + v > m, the last column takes on the simple form 

, _ r i , ( o ~ \ I 1 < i < s + t 

^t+2,s+t+lCis+t-lA^2,rn) i{s+t+l) = \ * „ ft __w- , _ ^ . , n " ( 21 ) 




Therefore, developing the determinant of -4f_J_ 2 m with respect to the last column, 

H+2,m dIlu ^l.m 



we are able to relate the determinants of Affn „ and Af'Z as 



det A^ m = x 2it - p) *2 P+ „j det ^ . (22) 

We are left with the task to specify the behaviour of the function g with respect to 
the "reduction" of the first two variables 

f&m = * l - m+n+2 x l ~ m+n+2 °t 9ft ■ (23) 



Assembling the two factors (|22|) and (|23|) , we obtain, in terms of the parameteriza- 
tion (1131) 



fiWW i_ -s2s+T+l + I 2(s-p+l)+T-T'+/i~- I (OA\ 

^2s+2+r,2t+T> — \ l ) u 2s+r \ X u 2p+u I ^C2s+T,2t+r' • 

K p=0 



We are now in the position to compare our general construction ( p4|) with the 
recursive equation for the Q-polynomials of the SU(3) 2 -HSG model (g). We read 
off directly the following restrictions 

v = ( and C = t' - p, - 1 . (25) 
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A further constraint results from relativistic invariance, which implies that the 
overall power in all variables of the form factors has to be zero for a spinless 
operator. Introducing the short hand notation \F®\ for the total power, we have to 
evaluate 



Combining (|2f 



Q 

and 



2S+T,2t+T' 



92s+T,2t+r' 



+ 



det A.2 S+Tt 2t+T' 



with the explicit expressions [det A:. 



2s+T,2t+r'\ 



(26) 
s{2t + v) + 



= l{l-m + p)/2 + m(p-m)/2 and [Q^f T % t+T ,] = l(l-l)/2-m(m-l)/2, 
we find the additional constraints 

l+r-i> 



and 



TV = t'(v — 1) 



(27) 



Collecting now everything we conclude that different solutions to the form factor 
consistency equations can be characterized by a set of four distinct integers. Assum- 
ing that each solution corresponds to a local operator, there might be degeneracies 
of course, we can label the operators by /i, v, r, r', i.e. O — > Off'",, such that we can 
also write Q^i instead of O^f"* 1 * ■ Then each Q-polynomial takes on the general 
form 

n O\2s+T,2t+r' _ n O^,\2s+T,2t+r' _ n ^ v 
V — V ' — V2s+r ! 2f+r' 

and the integers /i, u, r, r' are restricted by 

rv + r'fM = tt 1 , 2 + n> r, 2 + v > r' 



t'~ 1— U,T+1— fM 
92s+T,2t+r' 



det^ +rj2t+T , (28) 



(29) 



We combined here (^5]) and ( p7|) to get the first relation in (|29|) . The inequalities 
result from the requirement in the proof which we needed to have the form ([21]) . 
We find 12 admissible solutions to (|2l^), i.e. potentially 12 different local operators, 
whose quantum numbers are presented in table 1. 
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1/10 
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1/10 
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-1/2 





1/10 
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-1/2 





1/10 
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-1/2 





1/10 
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-1/2 


1/10 


1 





2 








-1/2 


1/10 


1 





1 








-1/2 


1/10 
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-1/2 


-1/2 


1/10 
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-1/2 


-1 
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1 














-1/2 


* 



Table 1: Operator content of the S'C/(3)2-HSG model. 
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Comparing with our previous results, we have according to this notation F °> ' ' = 

F fi\2s,2t^ pO°'*\2s,2t+l _ p£\2s,2t+l &nd pO\ ;*|2«,2f+l ^ _p6|2s+2,2t+2_ The ^ twQ 

solutions are only formal in the sense that they solve the constraining equations 
(p9|), but the corresponding explicit expressions turn out to be zero. 

In summary, by taking the determinant of the matrix (|TT|) as the ansatz for the 
general building block of the form factors, we constructed systematically generic 
formulae for the n-particle form factors possibly related to 12 different operators. 



4 Momentum space cluster properties 

Cluster properties in space, i.e. the observation that far separated operators do 
not interact, are quite familiar in quantum field theories p3] for a long time. In 
1+1 dimensions a similar property has also been noted in momentum space. For 
the purely bosonic case this behaviour can be explained perturbatively by means 
of Weinberg's power counting theorem, see e.g. || This property has been 

analysed explicitly for several specific models [21], [T^, It states that whenever 



the first, say k, rapidities of an n-particle form factor are shifted to infinity, the n- 
particle form factor factorizes into a k and an (n — «)-particle form factor which 
are possibly related to different types of operators 

• • • A) ~ *f(0 Xl . . . , e K )F°: K (e K+1 , . . . a) . (so) 

For convenience we have introduced here the operator 

6 

A— oc 



T i» ■= & n n (si) 



which will allow for concise notations. It is composed of the translation operator 
T„ A which acts on a function of n variables as 



r a A /(0i,...A,...A) ^ f(6i,...,e a + \,...,e n ) . (32) 



Whilst Watson's equations and the residue equations, see e.g. f| [7|, |16 



arc 



operator independent features of form factors, the cluster property captures part of 
the operator nature of the theory. The cluster property ( |3~0D does not only constrain 



the solution, but eventually also serves as a construction principle in the sense that 



* There exists also a heuristic argument which provides some form of intuitive picture of this 
behaviour ]l5| by appealing to the ultraviolet conformal field theory. However, the argument is 
based on various assumptions, which need further clarification. For instance it remains to be 
proven rigorously that the particle creation operator V^(9) tends to a conformal Zamolodchikov 
operator for 9 — > oo and that the local field factorizes equally into two chiral fields in that situation. 
The restriction in there that lim^oo Sij(8) = 1, for i being a particle which has been shifted and 
j one which has not, excludes a huge class of interesting models, in particular the one at hand. 
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when given F® we may employ fl3"0|) and construct form factors related to O' and 
O" . Hence, (|30D constitutes a closed mathematical structure, which relates various 
solutions and whose abstract nature still needs to be unraveled. 

We shall now systematically investigate the cluster property (|30|) for the SU (3)2- 
HSG model. Choosing w.l.g. the upper signs for the particle types in equation (^|), 
we have four different options to shift the rapidities 

sr±\ TpO\lx+,mX — TT^ TpO\lx+,mX — /qq\ 

1 l,K<l r n — 2 K +l<l,n r n \°°) 

q-±X pO\lx+,mx- _ pO\lx+,mx- ('34") 

which a priori might all lead to different factorizations on the r.h.s. of equation 
(|30|). The equality signs in the equations fl33|) and ( |3~4D are a simple consequence 
of the relativistic invariance of form factors, i.e. we may shift all rapidities by the 
same amount, for O being a scalar operator. 

Considering now the ansatz (^) we may first carry out part of the analysis for 
the terms which are operator independent. Noting that 

and T± A F+;(0) ~ ( ^ , (35) 



we obtain for the choice of the upper signs for the particle types in the ansatz 



i<j l<i<j<K K<i<j<l+m 



K-l A«(l-i) 

<j k (xi,...,x k )^~ e 2 

°1-k{Xk + 1,-,Xi)2 

m — l + K \K(l — m—l) 

cr K (xi,...,x K ) 2 e 1 

K 

<T„- K (x K + \,...,X n )'Z 



n+l—K<m,n 



][[ 

i<j 



n i-'" "■((',,) n f~ 

l<i<j<n—K n—n<i<j<n 



K \K,(m — l — l) 

o- n _ K (xi,...,a: n - K )?e 2 

l — 771+ K 

<J K (x n+ i- K ,...,X n ) 2 

k A/t(l-m) 
ff m -«.(^l+lr"A-);)^e 2 
K— m 

cr re (x n+ l_ re ,...,a;n) 2 



The remaining cases can be obtained from the equalities ( |33| ) and ([34]). Turning 
now to the behaviour of the function g as defined in (|T2|) under these operations, 
we observe with help of the asymptotic behaviour of the elementary symmetric 
polynomials ( ffTD and (|72|) 



r ii<2 9Lm = [e ±AK a K (xi, . . . , x K )o-;_ K (a; K+ i, . . . , Xi)) : '2 ' e (<r m ) " 2 



I — m-^rfi 



(36) 



n+1— ft<m,n"i m 



,±Ak 



Cftl^n+l-K, • • • j %n)0'm-ti{Xl+l, ■ ■ ■ > ^ti-k)] 2 -(37) 



In a similar fashion we compute the behaviour of the determinants 



T A Hpt >4' 



l,m 
l.m 



T A Hpt A^ v 

1 n+l-2K-£<m,n ueL ^l,m 



^ (2K+e) ( ( r 2 ^) i ((-l) t a;)«+« 1 -^ det J3* 



2t+!^ 



det ^ 



Ik— £ ,m 



e As( 2K+5 ) (a 2K+? ) s det .A; 



7"-A Hpt AW — ((—■\Yrr + ) K+ ^ Hpt A^ v 



Z,m— 2re— £ 



(39) 
(40) 
(41) 
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We have to distinguish here between the odd and even case, which is the reason for 
the introduction of the integer £ taking on the values or 1. 




T, 



rr 

V 0,1/ 



X 



T. 



X 



a: 



X 




X 



o 

•\^o, i/ 



X 



X 



X 



T 



X 



Figure 1: Interrelation of various operators via clustering. In this figure we use the 
abbreviations T x = ^ 2 «+i<i> T 2 = T i,L+i<v T 3 = T n-2 K <m,n> T 4 = T~} 2K<mn . We also 
drop the 2s and 2t in the subscripts of the (9's. Xhe Tj on the links operate in both 
directions. 

Collecting now all the factors, we extract first the leading order behaviour in A 



2 l,K<l r 2s+T,2t+r' e 



-\k(±V+T 



r (l=Fl) ' 



- i n+l-K<m,n- 1 2s+r,2t+r' c 



(42) 

Notice that, if we require that all possible actions of T^ b x should lead to finite 
expressions on the r.h.s. of (^), we have to impose two further restrictions, namely 
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t' > v and r > fi. These restrictions would also exclude the last two solutions from 
table 1. We observe further that F 22 tends to zero under all possible shifts. Seeking 
now solutions for the set fi, u, r, r' of (|52|) which at least under some operations leads 
to finite results and in all remaining cases tends to zero, we end up precisely with 
the first 9 solutions in table 1. 

Concentrating now in more detail on these latter cases which behave like 0(1), 
we find from the previous equations the following cluster properties 

1 l,2K.+£<l r 2s+T,2t+r> r 2K+£,0 r 2s+T-2K-£,2t+r' \^°) 

rr-\ t?h,v p0,0 pfi,u (AA\ 

2 l,2re+£<r 2s+t ,2t+u r 2K+£,0 r 2s+r-2K~£,,2t+u K**) 



2 n+l-2K-£,<m,,n r 2s+T,2t+T' r 2s+t ,2t+r' -2n-$, r 0,2k+£ J 

1 A ~j P 0,0 ( AP,\ 

2 n+l-2K-^<m,n r 2s+fj,,2t+T / r 2s+^,2t+r' -2n-^ r 0,2k+£ • \ w ) 



We may now use (^)-(|4q) as a means of constructing new solutions, i.e. we can 
start with one solution and use ([13|)- (|46|) in order to obtain new ones. Figure 1 
demonstrates that when knowing just one of the first nine operators in table 1 it is 
possible to (re)-construct all the others in this fashion. 

4.1 The energy momentum tensor 

As we observed from our previous discussion the solution F 22 is rather special. In 
fact this solution is part of the expression which in [l6j was identified as the trace 
of the energy momentum tensor 

Q e| 2s+ 2 M = ^(^3) e -(t + i) Vi(TiF M 22t+2 (47) 

The pre-factor <Jia\ will, however, alter the cluster property. The leading order 
behaviour reads now 

T ±A p e\2s,2t T ±A F B\2s,2t A(l-«/2) / 4 q\ 

2 l,K<2s r 2 n+l-K<2t,n r e ■ 

We observe that still in most cases the shifted expressions tend to zero, unless 
k — 1 for which it tends to infinity as a consequence of the introduction of the 
o\0\. There is now also the interesting case k — 2, for which the A-dependence 
drops out completely. Considering this case in more detail we find 



j-±\pQ\2s,2t ^ F e|2,O F 0|2s-2,2t x I cr 1 {x 3 ,...,X2B+2t) /^q\ 



rj-±X pe\2s,2t ^ jpO\2s,2t-2-pG\0,2 x I cr 1 (x 1 ,...,x 2s +2t-2) /j-q\ 

n+l-K,n | ai(xi,...,X2a) 

o\(x-i,...,X2s+2t-2) 



Vl(x 2s +1 


...,X2s+2t) 


<7l(x3,.. 


,X2s+2t) 


0-\{x 2a + l 


...,X2s+2t) 


CTl(x3,.. 


,X2s+2t) 


o-i (xi, 


■ ■,X2s) 


(J\(xi,...,X2s+2t-2) 


<Tl(xi, 


■;X2s) 
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Note that unless s = 1 in (^9) or t = 1 in (|50| ) the form factors do not "purely" 
factorize into known form factors, but in all parity breaking factor emerges. 

We now turn to the = 2s or k = 2t for which we derive 

r ±x F ei2s,2t „ r ±x_^ F e\2 S ,2t „ e x { 2-t- s) (51) 

We observe that once again in most cases these expressions tend to zero. However, 
we also encounter several situations in which the A-dependence drops out altogether. 
It may happen whenever t = 2, s = or s = 2, t = 0, which simply expresses the 
relativistic invariance of the form factor. The other interesting situation occurs for 
t — 1, s — 1. Choosing temporarily (in general we assume m_ = m + ) if®' 0,2 = 
2Txm 2 , m = m_ = m+e 2G//7r , we derive in this case 



T A F e|2,2 _ ^2 ^2 /c 9 \ 

In general when shifting the first 2s or last 2t rapidities we find the following 
factorization 

<-r±\ rp@\2s,2t q-±X r?e\2s,2t p9|2s,0 pe|0,2t fco\ 

1 l,2s r 2 n+l-2t,n r r r • \ oc> ) 

This equation holds true when keeping in mind that the r.h.s. of this equation 
vanishes once it involves a form factor with more than two particles. Note that 
only in these two cases the form factors factorize "purely" into two form factors 



without the additional parity breaking factors as in (49) and 



5 Identifying the operator content 

Having solved Watson's and the residue equations one has still little information 
about the precise nature of the operator corresponding to a particular solution. 
There exist, however, various non-perturbative (in the standard coupling constant 
sense) arguments which provide this additional information and which we now wish 
to exploit for the model at hand. Basically all these arguments rely on the as- 
sumption that the superselection sectors of the underlying conformal field theory 
remain separated after a mass scale has been introduced. We will therefore first 
have a brief look at the operator content of the Gfe/6 r (l) r -WZNW coset models 
and attempt thereafter to match them with the solutions of the form factor consis- 
tency equations. For these theories the different conformal dimensions in one model 
can be parameterized by two quantities | 25| : a highest dominant weight A of level 



smaller or equal to k and their corresponding lower weights A obtained in the usual 
way by subtracting multiples of simple roots aj from A until the lowest weight is 
reached 
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Here h is the Coxeter number of G and p the Weyl vector, i.e. the sum over 
all fundamental weights. Denoting the highest root of G by if}, the conformal 
dimension related to the adjoint representation A(^,0) is of special interest since 
it corresponds to the one of the perturbing operator which leads to the massive 
HSG-models. Taking the length of ip to be 2 and recalling the well known fact 
that the height of ip, that is ht(ip), is the Coxeter number minus one, such that 
(ip ■ p) — ht(-ip) = h — 1, it follows that A ^'°^ is a unique operator with conformal 
dimension A(ip,0) = h/(k + h). Note that uniqueness demands in addition that 
we do not take the multiplicities of the A-states into account. For SU(3)2 the 
expression ([54]) is easily computed and since we could not find the explicit values 
in the literature we report them for reference in table 2. 



A\A 


Ai 


A 2 


Ai + A 2 


2Ai 


2A 2 


A 


1/10 


1/10 


1/10 








A — ai 


1/10 


* 


1/10 


1/2 


* 


A - a 2 


* 


1/10 


1/10 


* 


1/2 


A — (X\ — Ct 2 


1/10 


1/10 


3/5 


1/2 


1/2 


A - 2«! 


* 


* 


* 





* 


A-2a 2 


* 


* 


* 


* 





A — 1a\ — a 2 


* 


* 


1/10 


1/2 


* 


A — a\ — 2«2 


* 


* 


1/10 


* 


1/2 


A - 2a x - 2c*2 


* 


* 


1/10 









Table 2: Conformal dimensions for A ( A,A ) in the 5C/(3) 2 /C/(l) 2 -coset model. 



Turning now to the massive theory, a crude constraint which gives a first glimpse 
at possible solutions to the form factor consistency equations is provided by the 
bound [0 

~F^"^(9 h ...,e n )].<£P . (55) 

We introduced here lim^.^oo . . . , 6 n ) =: const exp([/(6 l i, . . . , n ))i0i) as abbre- 
viation and denote the conformal dimension of the operator O in the ultraviolet 
conformal limit by A . We use the notation [ ]± when we take the limit in the 
variable X{ related to the particle species p i = " ± " , respectively. For the different 
solutions we constructed, we report the asymptotic behaviour in table 1. When 
we are in a position in which we already anticipate the conformal dimensions the 
bound (p5|) will severely restrict the possible inclusion of factors like a\, a\, a{ , cr^, 
which as we mentioned above may always be added since they trivially satisfy the 
consistency equations. 

More concrete and definite values for A° are obtainable when we exploit the 
knowledge about the underlying conformal field theory more deeply. Considering 
an operator which in the conformal limit corresponds to a primary field we can of 
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course compute the conformal dimension by appealing to the ultraviolet limit of 
the two-point correlation function 

(0^)0,(0)) =E^r 2A - 2A '- 2 ^ (O fc (0)> + . . . (56) 
k 

The three-point couplings are independent of r. In particular when assuming 
that is the smallest conformal dimension occurring in the model (which is the 
case for unitary models), we have 

l r \ 1/2A °" 

]im(O(r)O(0))~r-^° forr« UT) ) • (57) 

\C A o A o A o" {U )) 

Here O" is the operator with the second smallest dimension for which the vacuum 
expectation value is non-vanishing. Using a Lorentz transformation to shift the 
0(r) to the origin and expanding the correlation function in terms of form factors 
in the usual fashion 

oo °9 °9 in in / n \ 

(O(r)O'(O)) = E E S -^T^^ -rE^cosh^ 

x F^-^(9 1 ,...,9 n ) (Ff IM-^,... A))*- (58) 

we can compute the l.h.s. of (|5T| ) and extract A° thereafter. The disadvantage to 
proceed in this way is many-fold. First we need to compute the multidimensional 
integrals in (|58| ) for each value of r, which means to produce a proper curve requires 
a lot of computational (at present computer) time. Second we need already a 
relatively good guess for A°. Third for very small r the n-th term within the sum 
is proportional to (log(r)) n such that we have to include more and more terms in 
that region and fourth we need the precise values of the lowest non- vanishing form 
factors, i.e. in general vacuum expectation values or one particle form factors to 
compute the r.h.s. of (0). However, the lowest non- vanishing form factor can be 
of an arbitrary particle number and one may still extract the value of A . 

A short remark is also due concerning solutions related to different sets of /x's. 
The sum over the particle types simplifies considerably when taking into account 
that form factors corresponding to two sets, which differ only by a permutation, lead 
to the same contribution in the sum. This follows simply by using one of Watson's 
equations || 0, || ||], which states that when two particles are interchanged we 
will pick up the related two particle scattering matrix as a factor. Noting that the 
scattering matrix is a phase, the expression remains unchanged. 

5.1 A-sum rules 

Most of the disadvantages, which emerge when using ([57]) to compute the con- 
formal dimensions, can be circumvented by formulating sum rules in which the 
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r-dependence has been eliminated. Such type of rule has for instance been formu- 
lated by Smirnov [21] already more than a decade ago. However, the rule stated 



there is slightly cumbersome in its evaluation and we will therefore resort to one 
found more recently by Delfino, Simonetti and Cardy |15 . In close analogy to the 



spirit and derivation of the c-theorem [F25J these authors derived an expression for 
the difference between the ultraviolet and infrared conformal dimension of a primary 
field O 

oo 

A° uv - A° ir = - — J r (e(r)O(O)) dr . (59) 

Using the expansion of the correlation function in terms of form factors (|58D we 
may carry out the r-integration in (1591) and obtain 



1 oo ™ ™ MM 

a<L-a° = -^E E f f x " 71 



x 



2 (O) n=lMl ... Mn _{ o _4 n \(2ity (e?=i ™ H cosh^) 
F?^"»»(9 1 ,...,9 n ) (F^-^(e u ...,e n )Y . (60) 



Notice also that unlike in the evaluation of the c-theorem, which deals with a 
monotonically increasing series, due to the fact that it only involves absolute values 
of form factors, the series (|59D can in principle be alternating. Before the concrete 
evaluation of the expression ( pD| ) for the various solutions we constructed for the 
SU(3) 2 -HSG model, we should pause for a while and appreciate the advantages 
of this formula in comparison with (|57D. First of all, since the r-dependence has 
been integrated out we only have to evaluate the multidimensional integrals once. 
Second the evaluation of (|60"D does not involve any anticipation of the value of A° . 
Third one is very often in the comfortable position that despite the fact that the 
vacuum expectation value occurs explicitly in the sum rule, its explicit form is not 
needed. Once it is non- vanishing, the next lowest non-vanishing form factor may be 
normalized such that (O) cancels from the whole expression. For a singular vacuum 
expectation value the value of the remaining integral guarantees that the sum rule 



maintains its form as for instance discussed in [27|. Finally and most important 
fourth, the difficulty to identify the suitable region in r which is governed by the 
(logr) n behaviour of the n-th term in the sum in fl5"8|) and the upper bound in (|56| ) 
has completely disappeared. 

There are however little drawbacks for theories with internal symmetries and for 
the case when the lowest non-vanishing form factor of the operator we are interested 
in is not the vacuum expectation value. The first problem arises due to the fact 
that the sum rule is only applicable for primary fields O whose two-point correlation 
function with the energy momentum tensor is non-vanishing. Since in our model 
the n-particle form factors related to the energy momentum tensor are only non- 
vanishing for even particle numbers, we may only use it for the operators Oq^, 
0\ o and 2 2) where the latter operator is plagued be the second problem. 
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We will now compute the sum rule for the operators Oq$, O ' 2 , O^'q up to the 
6-particle contribution. We commence with the two particle contribution which is 
always evaluated effortlessly. Noting that 

F®(6) = -2™ 2 sinh(#/2) (61) 

and the fact that A^ r is zero in a purely massive model, the two particle contribution 
acquires a particular simple form 

oo 

(A°) (2) = -r±=r f (F? ++ (28)Y . (62) 



4tt(C») J coshfl 



Using now the explicit expressions for the two-particle form factors (74]), we imme- 
diately find 

( A o °;S)(2) = ( A og;i)(2) = ( A ^:o°)(2) = 1/8 (63) 

We recall, see also ||16|| , that in the limit a — > oo we obtain two copies of the 
thermally perturbed Ising model. This means that in the sum over the particle 
types in (|60|) there will be no contributions from terms involving different types of 



(Ato 


)(4) 


= 0.091 


^7 


(A°o.'o 


)(6) 


= 0.1004, 


(A^ 1 


)(4) 


= 0.08c 4 


iO 


(A°0°2 


)(6) 


= 0.0895, 


(A^:<» 


)(4) 


= 0.08< 


SO 


(A°2.0 


)(6) 


= 0.0895 



particles. We only obtain two equal contributions, namely 1/16 from F^ ++ and 
F% , such that the operator Oq'q plays the role of the disorder operator, as we 
expect. 

To distinguish the operators 0®^, O^l, G 2 o fro m each other we have to proceed 
to higher particle contributions. At present there exist no analytical arguments for 
this and we therefore resort to a brute force numerical computation. 
Denoting by (A )*™) the contribution up to the n-th particle form factor, our nu- 
merical Monte Carlo integration^ yields 

(64) 
(65) 
(66) 

We shall be content with the precision reached at this point, but we will have a 
look at the overall sign of the next contribution. From the explicit expressions of 
the 8-particle form factors we see that for Cqo the next contribution will reduce 
the value for A. For the other two operators we have several contributions with 
different signs, such that the overall value is not clear a priory. In this light, we 
conclude that the operators O ' , O®'^, C^'o & h possess conformal dimension 1/10 
in the ultraviolet limit. Unfortunately, the values for the latter two operators do 
not allow such a clear cut deduction as for the first one. Nonetheless, we base our 
statement on the knowledge of the operator content of the conformal field theory 
and confirm them also by elaborating directly on ([77|) and Q5%p. 

§We employed here the widely used numerical recipe routine VEGAS p8| . Typical standard 
deviations we achieve correspond to the order of the last digit we quote. 
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Figure 2: Rescaled correlation function G° '^(R) := (oJ^W^C )) P art ( a ) and ( G e) {8) (R) ■= 
(G(i?)0(O)) part (b) summed up to the eight particle contribution as a function of 
R = rm. 
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5.2 A from correlation functions 



First of all we do not presume anything about the conformal dimension of the 
operator O and multiply its two-point correlation function fl58f ) by r p with p being 
some arbitrary power. Once this combination behaves as a constant in the vicinity 
of r = we take this value as the first non-vanishing three-point coupling divided 
by the vacuum expectation value of O and p/4as its conformal dimension. This 
means even without knowing the vacuum expectation value we have a rational to 
fix p, but we can not determine the first term in (|56|). Figure 2a) exhibits this 
analysis for the operator Oq^ up to the 8-particle contribution and we conclude 
from there that its conformal dimension is 1/10. For the other operators the figures 
look qualitatively the same. 

The results of the same type of analysis for the energy momentum tensor is 
depicted in figure 2b), from which we deduce the conformal dimension 3/5. Recall- 
ing that the energy-momentum tensor is proportional J29| to the dimension of the 
perturbing field this is precisely what we expected to find. 

Furthermore, we observe that the relevant interval for r differs by two orders of 
magnitude, which by taking the upper bound for the validity of ([)7]) into account 
should amount to Cj_ _L n Cz 3 j_/(C3 3 n Cj_ j_ j_) ~ O(10~ 2 ). Since to our knowledge 

!0 10 ,5 5 10 5 5 u 10 10 10 ' . . , 

these quantities have not been computed from the conformal side, this inequality 
can not be double checked at this stage. 



12/5 (n) 

R g 




0,10 



Figure 3: Rescaled individual n-particle contribution (R) to the correlation function. 
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In figure 3 we also exhibit the individual n-particle contributions. Excluding the 
two particle contribution, these data also confirm the proportionality of the n-th 
term to (log(r)) n . 

We have carried out similar analysis for the other solutions we have constructed 
and report our findings in table 1. We observe that the combination of the vacuum 
expectation value times the three-point coupling for these operators differ, which is 
the prerequisite for unraveling the degeneracy. 



6 Conclusions 

With regard to the main conceptual question addressed in this paper, we draw the 
overall conclusion that solutions of the form factor consistency equations can be 
identified with operators in the underlying ultraviolet conformal field theory. In this 
sense one can give meaning to the operator content of the integrable massive model. 
The quantity on which the identification is based is the conformal dimension of the 
operator. Naturally this implies that once the conformal field theory is degenerate 
in this quantity, as it is the case for the model we investigated, the identification 
can not be carried out in a one-to-one fashion and therefore the procedure has to 
be refined. In principle this would be possible by including the knowledge of the 
three-point coupling of the conformal field theory and the vacuum expectation value 
into the analysis. The former quantities are in principle accessible by working out 
explicitly the conformal fusion structure, whereas the computation of the latter still 
remains an open challenge. In fact what one would like to achieve ultimately is the 
identification of the conformal fusion structure within the massive models. 

It would be desirable to put further constraints on the solutions by means of 
other arguments, that is exploiting the symmetries of the model, formulating quan- 
tum equations of motion, possibly performing perturbation theory etc. 

Technically we have confirmed that the sum rule ( |59|) is clearly superior to the 
direct analysis of the correlation function. It would therefore be highly desirable 
to develop arguments which also apply for theories with internal symmetries and 
possibly to resolve the mentioned degeneracies in the conformal dimensions. 

It remains also an open question, whether the general solution procedure pre- 
sented in this manuscript can be generalized to the degree that the type of deter- 
minants presented will serve as generic building blocks of form factors. 

The specific conclusions for the SU(3)2-homogeneous Sine-Gordon model are 
as follows: We have provided a rigorous proof for the solutions of the form factor 



consistency equations which were previously stated in [flq| . In addition we found a 
huge number of new solutions. By means of the sum rule and a direct analysis of 
the correlation functions we identified the conformal dimension of these operators 
in the underlying conformal field theory. Considering the total number of operators 
present in the conformal field theory (see table 2) one still expects to find additional 
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solutions, in particular the identification of the fields possessing conformal dimen- 
sion 1/2 is outstanding. Nonetheless, concerning the physical picture presented for 
this model one can surely claim that it rests now on quite firm ground. After the 
central charge of the conformal field theory had been reproduced by means of the 
thermodynamic Bethe ansatz [TJ[] and the c-theorem in the context of the form 
factor program [[RJ , we have now also identified the dimension of various operators. 
In particular the dimension of the perturbing operator was identified to be 3/5. 
Acknowledgments: A.F. is grateful to the Deutsche Forschungsgemeinschaft 
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7 Appendix 

7.1 Elementary symmetric polynomials 

In this appendix we assemble several properties of elementary symmetric polyno- 
mials to which we wish to appeal from time to time. Most of them may be found 
either in or can be derived effortlessly. The elementary symmetric polynomials 
are defined as 

a k (xx, . . . ,x n ) = x h ...xi k . (67) 

h<...<i k 

They are generated by 



J] (x + x k ) = xn k <?k(xi, . . . , x n ) . (68) 

fe=l fc=0 

and as a consequence may also be represented in terms of an integral representation 

lr dz n 
a k ( Xl ,...,x n ) = — f ^-^Y[( z + Xk ), (69) 

Z.TXI J\z\=q Z k=1 

which is convenient for various applications. Here g is an arbitrary positive real 
number. 



With the help of ( j69|) we easily derive the identity 

a k (-x, x, xi, . . . , x n ) = a k (xi, ...,x n )- x 2 a k - 2 (xi, ...,x n ) , (70) 
which will be central for us. We will also require the asymptotic behaviours 



, , je* x a v (x 1 ,...,x 1J )<T k - ri (x v+1 ,...,x n ) forr]<k 

T ^ akiXl --- Xn) ~{e^ k (x u ...,x v ) for V >k (71) 
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and 



a k (x v+1 , . . . ,x n ) for rj<n-k 



T lv <J k (x U . . . , X n ) ~ < a k+v _ n ( xl ,..., Xrl )g n _ v ( Xrl+1 ,..., Xn ) f nrn ^ n L, ( ?2 ) 



which may be obtained from (|69|) as well 



7.2 Explicit form factor formulae 

Having constructed the general solutions in terms of the parameterization (|3]), it 
is simply a matter of collecting all the factors to get explicit formulae. For the 
concrete computation of the correlation function, it is convenient to have some of 
the evaluated expressions at hand in form of hyperbolic functions. 

7.2.1 One particle form factors 



F o 1>0 |+ = F o 0il | = F o 0il | = F o li0 |+ = fflfi = H0A (?3) 



7.2.2 Two particle form factors 

= z(0)tanhf, for O = O %°, 0$, O^ , (74) 

= H^e^F+r {d ), = (75) 

7.2.3 Three particle form factors 

P o|±±± _ H°> [ 1 1, :< - F^jOjj) o o o,i n i,o (7f .s 

7.2.4 Four particle form factors 

F e|++— -7rm 2 (2 + cosh(%)) 

~ 2cosh(^ 12 /2)cosh(^ 3 4/2) y min l " iiJ U ° 



^o%°l++- -(0o O ;o O )cosh(^/2 + ^/2) 

2cosh(^ 12 /2)cosh(^ 3 4/2) ^ " 



7.2.5 Five particle form factors 

0-0,1 r-r pA* 
^ lli<3 f m 

ni<i<i<5 cosh (^/ 2 ) 



f O| ±±±±± = U l<3 F ^ ( 9 f for 00,0^ O o,c, 5 ^1,^ O o,i (7Q) 
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7.2.6 Six particle form factors 



F e|++++— 7rm (3 + cosh(%)) 



o \°|++++~- {°0 0) ((^2 ) 2 + ^4 + ^2 /K + ^2 , , n ,„,„, , 

16cosh(0 56 /2)ni< i<J <4Cosh(%/2) ^ minK ' 



7.2.7 Seven particle form factors 

pO|±±±±±±± -H" '' n.;<j -^min J (^i) f n o,o n o,o n i,o n o,i / a9 N 



7.2.8 Eight particle form factors 

^61++ -7rm 2 (4 + J2i<j co8h(0 y )) cosh(fl 12 /2) ™ 



8 Il3<i<i<8 cosh(%/2) 



F e\++++ = 7rm 2 (q 4 ) 1/2 (a 1 at + )(4 + Ei<j cosh(%)) „ 

" 27 ni<l<i < 4 cosh(%/2) n 5 <Ki<8 cosh(%/2) min ^ v 

(84) 



P01++++++— -?rm 2 (4 + Ei<j cosh(%)) cosh(fl 78 /2) „ 

8 = an muo 7o\ ll^min \Pw (°5) 

8 lli< i<3 -< 6 cosJi(%/2J iKj 
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